Strangelets and the TeV-PeV cosmic-ray anisotropics 
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Several experiments (e.g., Milagro and IceCube) have reported the presence in the sky of regions 
with significant excess in the arrival direction distributions of Galactic cosmic rays in the TeV to 
PeV energy range. Here we study the possibility that these hotspots are a manifestation of the 
peculiar nature of these cosmic rays, and of the presence of molecular clouds near the sources. We 
propose that stable quark matter lumps or so-called strangelets can be emitted in the course of the 
transition of a neutron star to a more compact astrophysical object. A fraction of these massive 
particles would lose their charge by spallation or electron capture in molecular clouds located in the 
immediate neighborhood of their source, and propagate rectilinearly without decaying further, hence 
inducing anisotropies of the order of the cloud size. With reasonable astrophysical assumptions 
regarding the neutron star transition rate, strangelet injection and neutralization rates, we can 
reproduce successfully the observed hotspot characteristics and their distribution in the sky. 



Several experiments have reported strong anisotropy 
measurements in the arrival direction distributions of 
Galactic cosmic rays (CRs) in the TeV to PeV energy 
range (Super-Kamiokande, Tibet III, Milagro, ARGO- 
YBJ, and IceCube [TJ [2]). The data reveal the pres- 
ence of large scale anisotropies of amplitude ~ 0.1%. 
Smaller scale anisotropies of size ~ 10° — 30° are also 
detected with amplitude a factor of a few lower. Mila- 
gro has reported the detection at significance > 12a of 
two hotspots (regions with enhanced CR intensity) with 
amplitude ss 10 -4 , at a median energy of 1 TeV. ARGO- 
YBJ report similar excesses. IceCube observes localized 
regions of angular scale ~ 15° of excess and deficit in CR 
flux with significance ~ 5er around a median energy of 
20 TeV g]. 

The large scale anisotropy could be naturally explained 
by the diffusive transport of CRs within the Galactic 
magnetic fields 0]. On the other hand, the inter- 
mediate and small scale anisotropies are more difficult 
to explain. The main difficulty resides in the fact that 
the Larmor radius of particles in the TeV-PeV range 
is: r L « E/ZeB ~ 1.08 pc Z~ 1 (E/lPcV)(B/l fiG) -1 , 
where the magnetic field strength of the Galaxy is as- 
sumed to be B = lfiG (see [5] for a review). For par- 
ticles with tl -C l c , where l c = 10 — 100 pc is the co- 
herence length of the Galactic magnetic field (e.g., [5]), 
the propagation will be totally diffusive over a distance 
> l c . Neutrons would propagate rectilinearly, but their 
decay length around 10 TeV energies is less than 0.1 pc. 
These scales are far shorter than the distance of any close- 
by source capable of accelerating particles to PeV ener- 
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gies. Various phenomena, such as heliospheric modula- 
tion, neutron sources, nearby pulsars, peculiar structures 
of the local Galactic magnetic fields have been invoked, 
but none seem to give an obvious explanation [3]. 

In this work we study the possibility that the hotspots 
in the skymap are a manifestation of the peculiar nature 
of CRs, and of the presence of molecular clouds (MCs) 
near the sources. We propose that quark matter lumps 
or so-called strangelets could be produced and acceler- 
ated while a neutron star (NS) transitions to a quark 
star (QS). A fraction of these heavy particles would suf- 
fer spallation or electron capture in molecular clouds lo- 
cated in the immediate neighborhood of their source, 
and produce neutral fragments that would propagate 
rectilinearly without decaying further, hence inducing 
anisotropies of the order of the cloud size. 

STRANGELET PROPERTIES AND SOURCES 

Strangelets (also referred to as nuclearites) are sup- 
posed to be lumps of uds quark matter. According to 
the Witten hypothesis [5], ud matter is metastable and 
d-quarks decay by weak interaction, u+d — > u+s, to form 
more stable uds matter. These lumps could be formed in 
explosive events like a NS undergoing a phase transition 
to a QS as proposed in dark matter (DM)-driven scenar- 
ios [7] or in the high-density environments of a compact 
object merger event [5]. Direct searches are being con- 
ducted by, e.g., ground-experiments at the LHC (Alice 
and CMS experiment with the CASTOR calorimeters) 
or in space with the AMS-02 spectrometer. 

The mass number, A, of a stable strangelet can poten- 
tially range from a ~ 10 to A S> 10 10 . When considering 
a general non-zero strangeness-content, there is a nuclear 
stability valley for strangelets with baryonic number A 
and there is a poorly known minimum value of mass num- 
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ber A min w 10 — 600 [5] below which they are unbound. 
Typical values of strangelet binding energy are currently 
uncertain but supposed to be E/A ~ MeV — GeV en- 
ergies. There is not much information on their possi- 
ble charge value Z but it should be small and (most 
likely) positive for finite lumps [TO]. Several models of 
strangelets exist that lead to various Zj A dependencies. 
For example, for ordinary strangelets, Z = 0.3A 2 / 3 , while 
for CFL (color-flavour-locked) strangelets Z ~ 0.3A 1 / 3 
[TU] , Even smaller charge-to-mass ratios are allowed 

Z/A 10" 2 - 10~ 7 . Experiments such as CREAM 

and AMS-02 will have the ability to perform a direct 
measurement of the charge, and infer estimates of Z/A. 

If strangelets were responsible for the observed 
hotspots, they should produce detectable air-showers. 
This is possible if the kinetic energy per nucleon content, 
Kn, satisfies Kn — K to t/A > lGeV. Measurements 
indicate a total kinetic energy of particles in hotspots, 
if tot ~ E ~ TeV-PeV, which implies A < 10 2 - 10 4 . 

Neutron stars have been suggested as possible acceler- 
ators of strangelets [11 . Strangelets could be produced 
for instance in the course of a NS to QS transition [12] . In 
such events, a fraction / e j of the gravitational energy re- 
leased can be injected into the expelled outer crust, lead- 
ing to total kinetic energies E C j ~ 4 x 10 50 (/ O j/10 -3 ) erg 
for standard NS mass and radius 13J. The Lorentz 
factor of the ejected mass can be of order T ~ 
22(/ cj /lO- 3 )(12km/i?,)(M,/1.5A/ ) 2 (lO^ 5 M /M ej ), 
for NS mass M*, radius i?*, and ejected mass M ej [T5] , 
Particles of mass number A could then gain energies of 
order E acc ~ 21 (A/10 3 )(r/22) TeV, the typical energy 
observed in hotspots. 

Accelerated stranglets may experience energy losses by 
interacting with the radiation field close to the NS, and 
with the baryonic and radiative backgrounds of the su- 
pernova (SN) envelope. Refs. [13] concluded that there is 
room for the escape of accelerated particles. The discus- 
sion can be adapted to our case, except that strangelets 
are likely to have higher binding energies [10 , which will 
further help the escape. Besides, old NS may have a 
higher chance to undergo a transition |13L 116) , negligible 
radiative fields and no surrounding SN envelopes. 

INTERACTION WITH MOLECULAR CLOUDS 
AND HOTSPOTS CHARACTERISTICS 

Once they have escaped from the source, strangelets 
diffuse in the magnetized interstellar medium (ISM). 
The trajectory of strangelets should be totally diffu- 
sive, even when assuming a low charge. They can 
reach the Earth on a timescale At = d^/(2D) ~ 6 x 
10 5 Z 1 / 3 (d s /lkpc) 2 (£;/20TeV)- 1 / 3 yrs, where d s is the 
distance to the source and the diffusion coefficient is 
set to D(E) = 1.33 x 10 28 # kpc [£:/(3Z GeV)] 1 / 3 cm 2 s" 1 , 
with -f/kpc = -ff/(lkpc) the height of the Galactic halo 



[17] , The ionization and the spallation timescales in the 
ISM (of average density riisu = 0.5 cm~ 3 ) read respec- 
tively r ion - 7 x 10 12 Z- 2 (£:/20TeV)yrs, and r spa n ~ 
4 x 10 5 (A/10 3 )" 2 / 3 (n ISM /0.5 cm" 3 )" 1 yrs 10,, implying 
that spallation should affect particles only mildly during 
their flight from sources located within 1 kpc. 

Spallation could however play a prominent role if the 
source is born in or near a molecular cloud. Molecular 
clouds are the densest regions of the ISM, and consist 
mainly of molecular hydrogen. Their typical radius 
in the Galaxy is Rmc ~ 20 — 50 pc, and their density 
riMC ~ 10 2-6 cm~ 3 . In such regions, the spallation 
fraction can exceed unity, reaching r spa n = T csc /r spa ii ~ 
7.5Z 1 / 3 (i?Mc/25pc)(n M c/10 3 cm- 3 )(A/10 3 ) 2 / 3 , with 
r csc the diffusion time of strangelets in the cloud. 
The electron capture rate for strangelets in clouds 
with free electron density ~ rj e riMC (with rj e <C 1) 
should be a fraction of the ionization rate, of order 
n on ~ 10 _5 Z 7 / 3 ?7 e (n MC /10 3 cm" 3 ). As strangelets are 
predicted to be more bound than standard nuclei, these 
estimates can be viewed as upper limits for spallation. 
For electron capture, it is possible that the large size of 
strangelets dominates the effects of the charge, implying 
a scaling in ~ A 2 / 3 , and the rates quoted here can be 
viewed as a lower limit . 

A fraction of strangelets undergoing spallation or elec- 
tron capture (similar to that quoted for regular ions [14] ) 
may generate neutral secondaries. The work of [10] sug- 
gests that a tiny parameter space exists where spalla- 
tion could lead to bound neutral strangelets. Neutral 
strangelets can then propagate rectilinearly to the Earth 
and produce a hotspot in the sky of the angular size of the 
MC, 9 MC ~ 14° (i?Mc/25pc)(d M c/200pc)- 1 , with d MC 
the distance of the MC to the observer. Note that this 
corresponds roughly to the size of the observed hotspots. 

The MC can radiate a total energy in neutral 
strangelets of E^c = V E cj ^Mc/W d s-Mc)] 2 , with 
l(d s _Mc) = d^_ MC c/(2D), the effective distance trav- 
elled by a diffusing particle over the rectilinear distance 
4-MCi separating the center of the source to the cen- 
ter of the MC. The factor 77 is a free parameter that 
accounts for strangelet production rate at the source, 
and the low strangelet neutralization rate in the MC. 
This expression is only valid for d s -MC > Rmc- If the 

source is located at the center of a MC, all the produced 

(2) 

particles diffuse in the cloud, and -E^c ~ 7 lE C3 - The 
total energy of neutral strangelets radiated by the MC 
can then be expressed over the whole range of d s -MC a $ 

EmC = + 1 I E MC\~ 1 - 

The excess signal in a solid angle < around 
one source can be defined as the following signal- 
to-noise ratio: cr <Q = A r s ,<o/(A r i so ,<o) 1 / 2 , where 

N B <(i = iMC^( a ! ^Htt^s-mc^^ 1 indicates the num- 
ber of events expected in a solid angle < Q from a source 
and -/Vj SOi< ^ = EJi SOjST A(a,6) the corresponding number 
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of events expected for an isotropic background. For a 
MC located at coordinates (a, 6), at distance d,MC, an d 
separated by d s -MC from the source, the signal at energy 
E can then be estimated as: 



a(E) = 



V 

£3/2 



d 4 r 2 
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(1) 

where A(a, S) [in m 2 s sr] is the exposure of an exper- 
iment in the direction (a, S), Ji SOiSI (E) is the observed 
cosmic ray flux at energy E, per steradian, and At is 
the diffusion time for particles to travel over a distance 
min(i?MC7 2i?MC + <^mc)- For a source located inside the 
MC, the luminosity in neutral strangelets radiated by 
the MC at E = 20TeV is of order L M c = E MC /At ~ 
3.5 x 10 40 ?7Z- 1 / 3 (i?Mc/25pc)- 2 erg/s. 

Figure [TJ presents contours of the value of a (Eq. [TJ 
for strangelets with Z = 1, A = 10 3 at E = 20TeV, 
as a function of g?mc and g? s _mc, for a MC of radius 
-Rmc = 25 pc. For each set of distances (rfs-MCi^Mc), 
the signal a is calculated for a solid angle 51 corre- 
sponding to an angle in the sky of min(3°, #mc)- This 
takes into account the minimum smoothing angle of 
the anisotropy analysis conducted by Milagro and Ice- 
Cube (2.1°/cos((5) for Milagro [1] and - 3° for IceCube 
[2]). Cosmic ray measurements indicate Ji so , S r(20 TeV) ~ 
5 x 10 _17 eV _1 s _1 m~ 2 sr -1 , and we chose an exposure 
of A(a, S) = 10 13 m 2 s sr, roughly corresponding to the 
Milagro exposure at 20 TeV, over 7 years of operation. 

In our calculation, we set the efficiency factor to r] = 
5 x 10~ 8 . Efficiencies in the range of 10 -8 < r/ < 1CT 7 
lead to reasonable values in terms of a (as a cx 77), what- 
ever the relative location of the source and the MC, 
and the distance to the MC. From Eq. [TJ one can in- 
fer the strong dependency of a on the distance between 
the source and the MC: a cx d~_ MC , when the source 
is at the border of the MC. On the other hand, Fig. [TJ 
shows that the value of a is relatively constant as long 
as the source is at a relatively central position inside the 
MC. This range of r\ thus implies that only MC within 
1 — 2 kpc, and only sources located inside the MC can 
produce a significant hotspot (note also that local MCs 
are found beyond c?mc 70 pc). 

Such low values of this effective parameter rj leave room 
for combined uncertainties in possibly low strangelet in- 
jection at the source, strangelet acceleration, and neu- 
tralization efficiencies in the cloud. All these are largely 
unknown but it is expected that the fraction of strangelet 
ejected mass at the source should be at most ~ 1 — 10% 
of the mass difference in the transitioning NS and QS 
configurations |16j . The neutralization rates in the cloud 
should be a small fraction of r spa n and rj on cx n e for spal- 
lation and electron capture respectively. 

Magnetic fields in MCs are known to scale approx- 
imately with gas density as n 1 / 2 relative to the mean 
Galactic magnetic field. The diffusion coefficient scales 
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FIG. 1: Particle excess significance a (Eq. [TJ, as would be 
observed by Milagro with 7 years of data, at E — 20 TeV, 
as a function of the distance of the MC to the Earth, d-MC, 
and the distance between the source and the MC, d s _MC, 
for strangelets with Z = 1, A = 10 3 and a MC of radius 
Rmc = 25 pc, source of luminosity L M c = r]10 40 erg/s, and 
an efficiency factor r\ = 5 x 10~~ 8 . Color bar: value of a, black 
lines: specific numerical values of a as indicated. 
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cx _B _1 / 3 4 2 ^ 3 in the Kolmogorov diffusion 
regime. Taking into account these stronger fields would 
thus only result in an order of magnitude difference in <r, 
and the variations could be absorbed by the uncertainty 
in the efficiency factor rj. 

Whether the NS-QS transition actually gives birth to 
a pulsar is unknown. Hence strangelet sources will not 
necessarily be found at the position of an active source, 
and one reasonable assumption would be that they are 
distributed as old NS. In general these are everywhere, 
including in MC. Indeed, the number density of old NS 
in the Galaxy is of order - 10 -4 pc -3 [TBI , and if 10% 
of them do not get kicks at birth (i.e., remain in the 
MC), their mean separation is 20 pc. Molecular clouds 
and gamma ray pulsars have a similar distribution on the 
sky, while radio pulsars cover a much wider age range and 
have a broader distribution. The fraction of old neutron 
stars that are required to undergo NS-QS transition in 
order to account for the handful of observed hotspots is 
~ 10~ 4 . The diffusion time of strangelets in MCs, i.e., 
the time over which each source will be observable, is of 
order At; ~ Z 1 / 3 375 yrs. This implies a NS-QS transition 
rate of order 3 x 10~ 7 yr -1 . 
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COMPARISON WITH DATA AND SIGNATURES 

Most observed hotspots could be produced by MCs 
in the Gould Belt (a star forming region concentrating 
many MCs, that forms a ring at a distance from the Sun 
of ~ 0.7— 2 kpc), at the location where NS-QS transitions 
may have occurred. Interestingly, the Milagro hotspot la- 
belled "Region A" [T] lies in the direction of the Taurus 
Molecular Cloud, the nearest star formation region lo- 
cated at 140 pc, and that covers ~ 100 deg 2 in the sky 
[19] . "Region 1" of IceCube [2] is also in the direction of 
a remarkable MC: the Vela Molecular Ridge, located at 
0.7 — 2 kpc distance, of size ~ 15° in sky [15] , 

It is difficult to predict whether strangelets could pro- 
duce air-showers conspicuously different from those from 
ordinary cosmic rays, mainly because their cross-section 
is not known. Preliminary hadronic simulations with 
EPOS and CONEX show that the strangelet heavy mass 
(expected to produce shallow showers with low fluctua- 
tions) and its high binding energy (with the oposite ef- 
fect) could also compensate each other to produce or- 
dinary cosmic ray showers [20 . It leaves room for the 
possibility that the cosmic rays observed by [U [2] actu- 
ally be strangelets, that were not identified (as long as 
A ^5 10 2 — 10 4 , see earlier discussion). IceCube also re- 
ported composition measurements from 1 to 30 PeV, that 
are compatible with ordinary cosmic rays |21j . The muon 
rate for strangelets is expected to increase in oc A 0A , 
which should not have a noticeable impact for our values 
of A. This again is consistent with the IceCube reports. 

Finally, our model predicts point-like cosmic ray 
sources farther than a few kpc. In Fig. [TJ the cut-off 
of the signal at distances due ^ 1 — 2 kpc stems from 
the limited resolution angle and sensitivity of the instru- 
ments, which set the smoothing angle for the anisotropy 
search. Indeed, as the angular size of the MC in the 
sky diminishes, the signal is diluted inside one angular 
bin, and cannot be distinguished from the noise. With 
better resolution and sensitivity (with the High Altitude 
Water Cherenkov Observatory, HAWC, for example), the 
excess signal could remain high at larger distances, and 
point-like sources could be spotted. 

This scenario could also lead to multi-messenger 
signatures in secondary neutrinos or gamma-rays. These 
particles could be produced directly at the source 
when strangelets are generated, or when particles 
undergo spallation in the clouds. In the former case, 
the secondary signal should be point-like. In the latter 
case, the whole MC could be illuminated, and the 
expected gamma-ray flux should be comparable to the 
product of a multiplicity factor multiplied by the hot 
spot flux. The multiplicity factor is of the order of 
the spallation fraction in the MC or r spa u 10. The 
predicted flux over solid angle ft sr is N 7 (> ITeV) ~ 
5.10- 12 (cr/10)(r spal i/10)(fi/10^ 5 )s- 1 cm- 2 . This should 



be compared with the expected flux of cosmic ray 
induced gammas. 

We discussed the possibility that strangelets acceler- 
ated in nearby NS-QS transitions, and then becoming 
neutral by spallation or electron capture in molecu- 
lar clouds, could explain the small-scale anisotropics 
observed by several experiments at TeV-PeV energies. 
With reasonable astrophysical assumptions regarding 
NS-QS transition rates, particle injection and neutral- 
ization rates, we can reproduce successfully the observed 
hotspot characteristics and their distribution in the sky. 
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